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Abstract 

Mass parameters for the relative and neck motions in fusion reactions of symmetric systems 
^'^Zr+^'^Zr, ^^''Pd+^^'^Pd, and ^'^^Ba+^^^Ba are studied by means of a microscopic transport modeh 
The shape of the nuclear system is determined by an equi-density surface obtained from the density 
distribution of the system. The relative and neck motions are then studied and the mass parameters 
for these two motions are deduced. The mass parameter for the relative motion is around the 
reduced mass when the reaction partners are at the separated configuration and increases with 
decrease of the distance between two reaction partners after the touching configuration. The mass 
parameter for the neck motion first decreases slightly up to the touching configuration and then 
increases with the neck width, and its magnitude is from less than tenth to several times more than 
the total mass of the system. The mass parameters obtained from the microscopic transport model 
are larger than the ones obtained from the hydrodynamic model and smaller than those obtained 
from the linear response function theory. The mass parameters for both motions depend on the 
reaction systems, but the one for the relative motion depends on the incompressibility of the EoS 
more obviously than that for neck motion. 
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I. INTRODUCTION 



In the study of a large amplitude collective motion in nuclear systems including heavy 
ion fusion reactions( especially the synthesis of superheavy nuclei) and fission, the macro- 
scopic model plays a very important role, in which a few collective degrees of freedom are 
usually used to describe the complex dynamic process. The potential energy surface, the 
mass parameter and the viscosity are the most important quantities in the macroscopically 
description of the large amplitude collective motion. Many works have been done for the 
study of these three quantities. The importance of the potential energy surface in the nu- 
clear large amplitude collective motion such as in fission and heavy ion fusion process is well 
known and up to now the calculation of potential energy surface based on the macroscopic 
and microscopic method seems to be quite successful [1, 2, 3, 4]. But the problem of kinetic 
energy ,i.e. the mass parameter, is far less settled than that of potential energies. As a 
matter of fact, the dynamical behavior of nuclear systems with large- amplitude collective 
motions, such as fusion and fission not only depends on the potential energy surface but 
also on mass parameter. Take a most simple example, let us consider the fusion process 
as a barrier penetration problem. In the WKB formula for barrier penetration the inertia 
tensor M(R) appears in the action integral S(L)=^| J^{2M{R){U{R) - E)Y''^dR along the 
fusion path L. The U(R) and E are the potential energy surface and the center of mass 
energy, respectively. The sign '-' and '+' refer to above and below the barrier energies. For 
simplicity only the internuclear distance R is taken to be the collective coordinate. The 
probability for penetrating the barrier(may be a multi-dimension barrier) is calculated as 
usual by P(L)=exp(-S(L)). In many studies the M(R) was simply taken to be the reduced 
mass. This simple approximation is only partly valid for fusion reaction of light nuclei. 
For fusion reactions of heavy nuclei, it is well known that not only the relative collective 
motion but also the neck motion becomes very important. The process for penetrating the 
barrier becomes more complicated and the M(R) can not be simply taken to be the reduced 
mass. Thus, the impact of the mass parameter on fusion process for heavy systems becomes 
incontestable. Recently renewed interest in the study of mass parameter is motivated by 
the study of the synthesis of superheavy nuclei[5, 6, 7, 8]. The qualitative and quantitative 
study of the mass parameter becomes more important. In the dynamic study of synthesis of 
superheavy nuclei the mass parameters for collective motions were studied by the cranking 
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model( or by the linear response function theory(LRFT) based on the cranking model) and 
the hydrodynamical model. In [5, 6], a large mass parameter for neck motion was found 
by means of the LRFT method leading to a restriction for a growing neck in dinuclear sys- 
tem and melting of the dinuclear system along the internuclear distance into a compound 
system. While in [7] a hydrodynamical mass parameter was adopted, then a subsequent 
shape evolution of pear shaped mono-nucleus formed with two incident ions could happen. 
Thus, the rather different mass parameters obtained by cranking and hydrodynamical mod- 
els lead to complete different description of the fusion mechanism for heavy systems due to 
different fusion paths. The cranking model was developed to calculate the mass parameter 
of nuclear systems microscopically[9, 10, 11, 12, 13, 14]. But the abrupt change of mass 
parameters near the level crossing makes difficulties for practical usage[15]. Whereas for 
the hydrodynamical model, the assumption of the irrotational flow of an incompressible and 
non-viscous fluid [7, 16, 17, 18, 19] is too simple for many practical purposes. Furthermore, 
in both approaches, one usually adopts a way that only one collective degree of freedom 
is taken to be a variable and the others are fixed in the calculation of mass parameters. 
The condition of the validity of this (static) treatment needs to be tested because in a real 
heavy ion fusion reaction all degrees of freedom of collective motion change dynamically and 
self-consistently. Therefore, it seems to be necessary to investigate the mass parameters for 
the collective motion in fusion reactions by a more realistic dynamic model. In this work 
for the first time we try to employ the microscopic transport model, namely the Improved 
Quantum Molecular Dynamics (ImQMD) model to investigate the mass parameters for col- 
lective motions in heavy ion fusion reactions and possibly to find out the conditions for the 
validity of the widely used static treatment in the calculations of mass parameter. 

The quantum molecular dynamics (QMD) model being successfully used in intermediate 
energy heavy-ion collisions was extended to apply to heavy ion collisions at energies near 
barrier by making a series of improvements [20, 21]. The main improvements introduced are: 
introducing the surface and surface symmetry energy terms in the potential energy density 
functional, introducing the phase space occupation constraint [22] which can be considered as 
an approximate treatment of anti-symmetrization, and a system size dependent wave packet 
width. Furthermore, considering the especially importance of a proper initial condition for 
the application of the microscopic transport model to low energy heavy ion reactions, in 
the ImQMD model the requirement of the projectile and target at initial time being at 
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ground states is enforced. And consequently the projectile and target can be stable for 
a long enough time and the structure effect in entrance can be taken into account partly. 
Thus, in principle, the dissipation, diffusion and correlation effects are all included without 
introducing any freely adjusting parameter. With the ImQMD model the fusion dynamics at 
energies near and above the barrier has extensively been studied. Firstly, the stable initial 
nuclei with good ground-state properties such as the binding energies, root-mean-square 
radii and density distribution for a series of nuclei from ^^O to ^ospi^ and 238 ^ ^qj-q realized 
within several thousands fm/c [20, 23]. The Coulomb barriers for many reaction partners 
can be described well. The model were applied to calculate fusion (capture) excitation 
functions for a series of fusion reactions including neutron-rich projectile and target [21]. 
We also extend the ImQMD model to study the capture cross sections for the reactions 
^^Ca+^^'^Pb,^^^!! leading to super-heavy nuclei [24]. The calculated results are in good 
agreement with experimental data. Thus, it is more confident to apply the ImQMD model to 
study the mass parameters for collective motions in heavy ion fusion reactions. In this work, 
the mass parameters for symmetric reactions of ^°Zr+^'^Zr, ^^°Pd+^^''Pd, and ^^'^Ba+^^^Ba 
are investigated. The paper includes following parts: In section II we briefly introduce the 
ImQMD model. In section III, we present the results of our study on the mass parameters for 
^°Zr+^'^Zr, ^^'^Pd+^^'^Pd, and ^^^Ba+^^^Ba. Finally, we give brief summary and discussion 
in section IV. 



II. BRIEF INTRODUCTION OF THE IMQMD MODEL 

In the ImQMD model as the same as the original QMD model, each nucleon is represented 
by a coherent state of a Gaussian wave packet [25, 26], 

where, rj, pj, are the center of i-th wave packet in the coordinate and momentum space, 
respectively, represents the spatial spread of the wave packet. Through a Wigner trans- 
formation, the one-body phase space distribution function for N-distinguishable particles is 
given by: 

/(^.P)-E7i-p|-i^-f(p-p.n, (2) 
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The density and momentum distribution functions of a system read 

p(r)= f /(r,p)dp = 5^p,(r), (3) 

g{p)= j /(r,p)rfr = ^^7,(p), (4) 

respectively, where the sum runs over all particles in the system. pj(r) and 5'i(p) are the 
density and momentum distributions of nucleon i: 

AW = ^^;^-p[-i^:^], (5) 

^^(P)=(^;^-P[ (6) 

where a^. and (Xp are the widths of wave packets in coordinate and momentum space, respec- 
tively, and they satisfy the minimum uncertainty relation: 

a,ap = |. (7) 

The propagation of nucleons under the self-consistently generated mean field is governed 
by Hamiltonian equations of motion: 

dH . dH 



Hamiltonian H consists of the kinetic energy and effective interaction potential energy. 



(8) 



H = T + Uioc + Ucoui. (9) 

Here, Ucoui is the Coulomb energy, and Uioc = j Vioc[p{j')]dr . Vioclpij")] is the nuclear 
potential energy density functional, which is introduced according to Skyrme interaction 
energy density functional with spin-orbit term omitted [21, 27]. It reads 

VioAp] = f ^ + + f^(Vp)^ + ^{p' - ^s{Vpy)6' + g/-^. (10) 

2 po 7 + 1 Po 2po 2po Po' 

Here, 6 = {pn — Pp)/{pn + Pp) is the isospin asymmetry and the p, p„, pp are the nucleon, 
neutron, and proton density, respectively. By integrating Vioc, we obtain the local interaction 
potential energy: 
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E E f'S + 9. E ( E f ) + Y E E (1 - . (") 



where 



fr,- - r,;^2 



/s^i = ^ - (^^)'' (13) 

and tj=l and -1 for proton and neutron, respectively. 

The Coulomb energy can be written as the sum of the direct and the exchange contribu- 
tion 

Ucoui = \j j pp(r)^^p,(r')rfrrfr' - (^) ' ^ / pfdv. (14) 

The collision term and the constraint for single particle occupation number are considered 
as the same as that in [20, 21, 22]. 

To have a proper initial condition ( with good properties of projectile and target nuclei) is 
of crucial importance for studying low energy heavy ion reactions by means of the transport 
model description. In this work we pay special attention to the initial condition. The 
procedure for making initial nuclei is similar as that in ref . [20] . The initial nuclei applied in 
the study of reaction process have good ground state properties such as binding energies and 
rms radii and furthermore, their time evolution is very stable remaining almost unchanged 
for a long enough time(in this work, it is taken to be 2000fm/c). At the same time they are 
required to be without spurious particle emission. For the self-consistency, we adopt the 
same effective nuclear potential energy density in making initial nuclei and the simulation 
of the reaction process when solving the Hamiltonian equation(8). Two sets of interaction 
force parameters used in the calculations and the corresponding properties of saturated 
nuclear matter are given in table I. In this work, the definition of a fusion event is taken to 
be the same as that in [20], in which a fusion event is defined operationally like in TDHF 
calculations that is that the fusion event is defined rather operationally as the event in which 
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the coalesced one-body density survives through one or more rotations of composite system 
or through several oscillations of its radius. 





a(MeV) /3(MeV) 
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w(MeV/m2) 


5r(MeV) 


7] c. 


,(MeV) 






3) K^{MeV) 


IQl 


-310 258 


7/6 


19.8 


9.5 


2/3 


32 


0.08 


0.165 


144 


IQ2 


-356 303 


7/6 


7.0 


12.5 


2/3 


32 


0.08 


0.165 


175 



TABLE I: the model parameters 



III. THE MASS PARAMETERS FOR COLLECTIVE MOTIONS IN HEAVY ION 
FUSION REACTIONS 

In the macroscopic model, the low energy heavy ion fusion reactions are generally de- 
scribed by five collective coordinates, namely the distance between the centers of mass of 
two colliding nuclei, the mass asymmetry degree of freedom, the neck parameter and defor- 
mation coordinates of two ends of the colliding system. Since in the early stage of reactions 
the mass parameters for relative collective motion and neck motion are weakly dependent 
on the deformations of two end parts[6], one often considers only three collective modes, i.e., 
the relative motion of two nuclei, the formation and rupture of neck and nucleon transfer 
between projectile and target. The mass parameters for those collective motion modes are 
expressed by a matrix and the corresponding collective kinetic energy is given by 

T = \Y.^^M)m- (15) 

Here, Mjj is the mass parameter matrix and qi is the i-th collective coordinate. The nuclear 
mass parameters refiect the properties of a nuclear system related to the collective kinetic 
motion, which depend on the time dependent change of the shape of nuclear system and may 
depend on the effective nucleon-nucleon interaction. For the symmetric reactions considered 
in this work the mass asymmetric motion plays a negligible role and thus we can only 
consider the relative motion of projectile and target and the neck motion in the following 
study. 



7 



A. The mass parameter for the relative motion of projectile(like) and target(like) 



Let us first investigate the mass parameter for relative motion of projectile and target 
nuclei. In order to study the collective motion with the microscopic transport model one has 
to first define the collective coordinate by means of the microscopic quantities. The distance 
between the centers of mass of two nuclei R can be expressed as 



R = A^ z ■ sign(z)f{f,p)dfdp. 



(16) 



-1, z>0; 



Here, sign{z) = < and the f{r,f)) is the distribution function in phase space. 

-1, z<0. 

The distance between the centers of mass of two nuclei R and the neck width A (which will 
be discussed in the next subsection) for an axial symmetric nuclear system are sketched in 
Fig.l. In addition, the system length L is also shown in Fig.l. Within the ImQMD model 
approach, the conjugate momentum of the relative motion of two nuclei is expressed as 



Pz ■ sign{z) f {f, pjdfdp. 



(17) 
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FIG. 1: (Color online) The schematic figure of deformation parameters of a reaction system 

The velocity of the relative motion of two nuclei 
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At 



(18) 



can be obtained numerically. The mass parameter for relative collection motion of two nuclei 
is then calculated by 

Mrr = ^. (19) 
R 
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This formula seems to be questionable when the R approaches to zero, but , in fact, the Pr 
also approaches to zero as the i? ^ 0. This situation only happens at the fusion reaction 
closing to the end, which is a not interesting case. 

In this work the fusion systems of ^°Zr+^°Zr, ^^'^Pd+^^°Pd, ^^^Ba+^^^Ba are chosen for 
studying mass parameters. They are all along the /5-stability line. The head on collisions at 
the incident energy of 1.1 times the height of Bass barrier [29] are considered for all systems. 
The incident energy is selected in such a way that the fusion process can be proceeded, 
whereas the compress or expand of nuclear matter does not appear. 
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FIG. 2: (Color online)The mass parameters for the relative collection motion as a function of R/Rq 
for systems of ^^Zr+^OZr, ^lOpd+^^^Pd and ^^SBa+i^SBa. 



Fig. 2 shows the dependence of the mass parameter Mrr in the unit of total mass of the 
reaction system Mq on the scaled distance between the centers of mass of two nuclei R/Rq 
for systems of ^^Zr+^^Zr, ^^opd+^^opd and ^^SBa+i^SBa calculated with IQ2. Here Rq= 
tq^Ai + y42)^^^fm is the radius of the spherical compound nucleus formed by two nuclei. 
The ro is taken to be 1.16 fm. The corresponding configurations of the reaction system at 
different stages are also shown in the figure. The structures appeared in the Mrr in the 
figure, especially for ^^^Ba+^^^Ba is mainly due to not enough number of the fusion events. 
From the figure one sees that when R > I.SRq the mass parameter Mrh approaches to 0.25 
total mass of the system, i. e., the reduced mass, whereas in the case of R < I.SRq the 
Mrji increases with decrease of the distance R. Another important feature is that the mass 
parameter for relative motion of two nuclei depends on the size of the reaction systems. The 



heavier system has the steeper increasing slope of Mjir with decrease of R. 

The incident energy dependence of mass parameters is also interesting. Here we only 
study the energy dependence of the Mjui. Fig.3 shows the mass parameter for relative 
motion for ^^'^Pd+^^^Pd at E=1.0, 1.1 and 1.5V b, where refers to the height of the 
Bass barrier. The figure shows that the mass parameter for relative motion depends on the 
incident energy. For lower energy like E=1.0Vs, the mass parameter for relative motion 
changes slowly since the reaction process is slow and the reaction system has time to change 
their shapes gradually and consequently the mass parameter for relative motion increases 
gradually from reduced mass to larger than the reduced mass. As the incident energy 
increases the reaction process becomes faster and there is less time to change the shape of 
the system and the mass parameter for relative collective can be close to the reduced mass 
for a period of time and then increases rapidly at very late time. So we see from the figure 
that at higher energy such as at E=1.5Vb the Mjui is roughly equal to reduced mass when 
R/R-o > 1-3. But we expect it could be the case only for not heavy nuclear systems and 
at very above the barrier energies. The time evolution of the mass parameters for relative 
motion at three energies is shown in the inserted figure of Fig.3. 
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R/RO 

FIG. 3: (Color online)The incident energy dependence of the Mrr for ^^"Pd+^^^Pd at 
E=1.0,1.1,1.5Vb- The inserted figure is the time evolution of the Mrr for the same system. 

Now let us investigate whether the mass parameter for relative motion depends on the 
nucleon-nucleon effective interaction. In Fig. 4 we make a comparison between the mass 
parameters Mrr calculated with two force parameter sets IQl and IQ2 (listed in Table I) 
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FIG. 4: (Color online)The comparison between mass parameters M/j/j calculated with force pa- 
rameter sets IQl and IQ2. 

for reaction systems of ^^Zr Zr and ^^^Ba +^^^ Ba. One sees that the M^jr calculated 
with IQl rises slower than that calculated with IQ2. We notice that the parameter set IQl 
gives relatively smaller incompressibility compared with the force IQ2 (see Table I). The 
difference between the Mnji calculated with IQl and IQ2 can be attributed to the difference 
in incompressibility of interaction parameters IQl and IQ2. The softer EoS needs less energy 
to make two nuclei even closer than touching configuration and thus the mass parameter for 
relative motion for the fusion reaction system decreases with the decrease of the stiffness of 
the EoS. 

In order to make comparison with other approaches, in Fig. 5 we show the results cal- 
culated with hydrodynamic model(Werner- Wheeler (W-W) method) [7, 18] and with the 
linear response function theory model(LRFT)[6], in addition to the microscopic transport 
model approach for the system of ^^'^Pd+^^'^Pd. The results obtained by the ImQMD model 
with IQ2 and by hydrodynamic model(W-W) are shown by solid and long dashed curves, 
respectively. The short dashed curve is the results by LRFT taken from [6] where the large 
fluctuation appears due to the level crossing. The general tendency of the changing 
with R/Ro calculated by the ImQMD model is similar with those obtained by hydrody- 
namic model and by LRFT. The magnitude of the mass parameter obtained by the ImQMD 
model is higher than that obtained by hydrodynamic model, but lower than that obtained 
by LRFT. 
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FIG. 5: (Color online)The comparison of mass parameter Mrr calculated by the ImMQD model 
with that by Werner- Wheeler method and the linear response function theory. 

B. The mass parameter for the neck motion in fusion reactions 

In order to study the mass parameter for the neck motion within microscopic transport 
model one has to first define a nuclear surface based on the density distribution of the 
reaction system. The shape of the nuclear system is then defined by the nuclear surface. 
The nuclear surface is usually defined by an equi-density surface of the reaction system. 
In this work this equi-density surface, namely the ps{z,t) for a axial symmetric system at 
time t is taken to be at the half normal density p = O.Spo- The neck width A is defined 
by 2ps{z = Zc,t) and the Zc denotes the position of the neck as shown in Fig.l. Thus, the 
distance of two nuclei and neck width are evolved with time simultaneously. Fig. 6 shows the 
correlation between the neck width A and distance between two nuclei R/ Rq. The contour 
plots of density distribution for ^^^Ba+^^^Ba at R/ Ro=1.9, 1.8, 1.6 are also shown in the 
figure, the lines in the contour plots correspond to p=0.5, 0.75, 1.0 po, respectively. One 
can find that the touching configuration is around R/ Ro=1.8, but not at R/ Ro=1.6 which 
is about the sum of the radii of two nuclei. This is because two nuclei are elongated due 
to the interaction between two nuclei. One can further find that the most fast increase of 
the neck width A happens around the touching configuration, for instance, in the system 
^^^Ba+^^^ Ba the A from 0.5fm increases to about 7fm when the R/ Rq decreases only from 
1.9 to 1.6. It means that the touching configuration is important for characterizing the neck 



12 



motion. 
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FIG. 6: (Color online)The correlation between neck width A and distance R/Rq for ^^^Ba+^'^^ Ba. 

Mass parameters, in fact, are related with the kinetic energies of collective motions. The 
collective kinetic energy for relative motion of two colliding nuclei and neck motion can be 
expressed by 



As is indicated in ref. [15] and [5] that for the mass symmetric reactions the mixing mass 
parameter involving the neck degree of freedom remains comparatively small. It means that 
the condition M^a/ V MhhMaa << 1 is approximately satisfied [5]. In order to test whether 



that the coupling between the velocities of the relative motion and neck motion is indeed not 
strong within the time period which we are interested in and thus, neglecting the crossing 
term in (20) is reasonable. Here we should mention that the velocities for neck and relative 
motion are strongly correlated when t>450fm/c seen from the figure. We should also notice 
that the neck is already well developed at that time and in this case it is difficult to well 



T — Tnn + Tra + ^AA — - 




(20) 



the crossing term in (20) can be eliminated or not, in Fig. 7 we show the time evolution of R 
and A for ^^^Ba+^^^Ba. It is seen from the figure that the A is a strongly peaked function 
of time, and the peak is at about 350 fm/c with very steep rising left side, while i? is a 
monotonically decreasing function of time. The R reduces considerably at 280 fm/c, which 
is much earlier than the time when the A reaches the peak value. From this figure we find 
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distinguish the relative and neck motion and the treatment of neglecting the crossing term 
adopted in this work only gives approximate results . 
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FIG. 7: (Color online)The time evolution of the A and R for ^^SBa+^^Sgj^^ r^^^ ^^j^^ 

curve is for 

A and the dashed curve for R. 

Now let us investigate the mass parameter for the neck motion in symmetric fusion 
systems. Instead of using the relation of Maa=Pa/A we calculate Maa through the relation 

Taa = ^MaaA^ (21) 

and 

Taa = Etot — Epot — Trr — Toth — E^xc- (22) 

Here Toth is the kinetic energy from the motion of the additional degrees of freedom except 
the relative motion and neck motion. Since in the early stage of fusion reaction the change 
of the deformations at two ends of colliding system is not strong, the energy corresponding 
to this motion could be neglected[6]. Thus the Toth is mainly from the crossing term of 
relative motion and neck motion. As is mentioned above this term is small, then the Toth 
should be small and could be neglected. The E(.xc is the internal excitation energy. As we 
know that the neck motion plays important role only near touching configuration when the 
internal excitation is weak, in addition the incident energy is selected to be low, the internal 
excitation energy should be much smaller than the collective motions. Thus, we neglect the 
Eexc in the present calculations for simplicity. The Epot is calculated by 

Epot{R) = Eu{R)-E,-E2. (23) 
14 




Here, £"12 (-R), Ei and E2 are the total energy of the whole system, the energies of the 
projectile(like) and target (like) part, respectively. The reactions ^°Zr+^°Zr and ^^^Ba+^^^Ba 
at incident energies equal to 1.08 times the height of the Bass barrier (Ecm=350 MeV) are 
selected to study the mass parameter for neck motion. Here the incident energy chosen is 
lower than that in the investigation of the mass parameter for relative collection motion 
in order to reduce the internal excitation energy. Fig. 8 (a) shows the time evolution of 
kinetic energies for relative and neck motions, and of the potential energy for ^'^'^Ba+^^^Ba 
at Ecm=350 MeV. The parameter set IQ2 is utilized. Here the initial time is taken at the 
time when the surface to surface distance of two nuclei equals to 22.5fm (correspond to 
R/-Ro=3.6 in Fig. 8(b)). One sees from the figure that the kinetic energy for relative motion 
decreases firstly and then approaches to zero with time increasing, while the kinetic energy 
for neck motion is about zero at beginning, then it increases and until 550 fm/c it saturates. 
The potential energy is zero at the infinite distance, then it increases as two nuclei approach 
with each other and it reaches a maximum value at about 300 fm/c, after then it reduces to a 
saturated value. In Fig. 8(b) we show the kinetic energies for relative motion and neck motion 
and the potential energy as a function of R/Rq. The energy for relative motion gradually 
decreases with R/Rq decreasing and finally equals to zero, while the energy for neck motion 
starts to increase quickly from zero to a saturated value from R/Rq ~ 1.9 to R/Ro ~ 1.3. 
We find that the magnitude of the kinetic energy for neck motion exceeds that of relative 
motion just after the touching configuration. The kinetic energies for both relative and neck 
motion are saturated at about R/Ro=1.3. It is also noticed that the potential energy also 
reaches a maximum value near the touching configuration. 

Having the kinetic energy for neck motion Taa and velocity A, we can calculate the mass 
parameter for neck motion through expression (21). Fig.9(a) presents the evolution of the 
mass parameter for neck motion with the relative distance between the centers of mass of 
two nuclei and Fig. 9(b) shows the mass parameters as a function of the neck width scaled 
by L, the total length of the colliding system, for the systems ^^Zr Zr, ^^^Ba +^^^ Ba, 
respectively. From Fig. 9 one sees that the Ma a reaches a small minimum near the touching 
configuration. Then it increases with decrease of R/Rq and increase of A. We notice the 
change in the slope of Maa with respect to R/Rq and A near the touching configuration, 
which is due to the sharp peak of A as shown in Fig. 7. The magnitude of the Maa is from 
less than tenth of the system mass {Ai + A2) to several times the system mass. Finally it 
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FIG. 8: (Color online) The kinetic energies for the relative and neck motions and the interaction 
potential energy for the reaction ^^^Ba+^'^^Ba at Scm=350 MeV (a) as a function of time and (b)as 
a function of R/iio • 

approaches to saturate values for both systems of ^^Zr Zr and ^"^^Ba +^^^ Ba. The figure 
shows that slopes for Maa vs A(R/Ro) for ^^Zr Zr and ^^^Ba +^^^ Ba are different 
and it implies the system dependence of the mass parameter for neck motion. We have also 
investigated the influence of the effective interaction on the mass parameter for neck motion 
by comparing the mass parameters calculated with IQl and IQ2 and we find that there is 
no big difference between them. It seems to us that the dependence of the mass parameter 
for neck motion on the effective interaction is weaker than that for the relative collection 
motion in nuclear fusion reactions along the /3-stability line. 
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FIG. 9: (Color online)The mass parameter Maa for the neck motion as a function of a)R/Ro and 
h)A/L for systems ^^^Ba +^38 ^nd ^^Zr+^^Zr calculated with IQ2. 

Finally, we make comparison between the mass parameter for neck motion obtained by 
means of the microscopic transport model and the hydrodynamic model as well as the 
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LRFT. Fig. 10 shows Maa/Mq as function of neck widths A at R/Ro around the touching 
configuration. Here, the sohd, dashed, and dotted hues denote the Maa calculated by 
the microscopic transport model, the hydrodynamic model and the results taken from ref. 
[6] calculated by LRFT, respectively. The general behavior of the Maa as a function of 
neck width is coincident for all results obtained by three different approaches and also is in 
agreement with other studies[5, 15], i.e. the Maa increases with A and finally it saturates 
when the neck is already very much developed. The results of this work are in between the 
results from hydrodynamic model and LRFT approaches. This is similar with the situation 
in mass parameter for relative motion as shown in Fig.5. 
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FIG. 10: (Color online)The comparison of the mass parameter Maa calculated by the ImQMD 
model with that calculated by hydrodynamic model and by LRFT taken from [6]. 



IV. SUMMARY AND DISCUSSION 

In summary, in this paper we employ the microscopic transport model, namely the Im- 
proved Quantum Molecular Dynamics model to study the mass parameters of collective mo- 
tions in fusion reactions for the first time. The head on reactions of ^°Zr+^'^Zr, ^'^^Pd+^^^Pd, 
and ^^^Ba-|-^^^Ba are selected. The essential difference of the present approach from other 
methods is that within a microscopic dynamical model approach, the shape of the reaction 
system is determined by the time dependent density distribution of the system, so that the 
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distance between the centers of mass of two nuclei and the neck width change with time 
self-consistently during the reaction process. It is found that the feature of the relative 
and the neck motion is rather different, and the A and the R are coupled weakly, as seen 
from Fig. 7. So that we can study the M^a and M^/j individually. Then we investigate 
the mass parameters for the relative and neck motion for three systems. We find that the 
mass parameter for relative motion between two nuclei approaches the reduced mass when 
the reaction system is at the separated configuration and after contact of two partners it 
increases with decrease of the distance between two centers of mass. The mass parameter 
for neck motion has a small minimum near touching configuration and after touching con- 
figuration it increases with the increase of neck width(or the decrease of relative distance 
between two nuclei). Its magnitude is from less than the tenth to more than several times 
the total mass of the system. The general tendency of the dependence of M^ji on the R/-R0 
and Maa on A is similar with those obtained by the hydrodynamic model and the LRFT. 
The magnitude of mass parameters obtained in this work is larger than the ones obtained by 
the hydrodynamic model and smaller than those obtained by LRFT. Both Mur and Maa 
depend on the reaction systems. And the influence of the effective interactions on the mass 
parameters is obvious for the mass parameters M^ji but not for Maa- 
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